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Abstract. Conductivity equation is studied in piecewise smooth plane do- 
mains and with measure- valued current patterns (Neumann boundary values). 
This allows one to extend the recently introduced concept of bisweep data to 
piecewise smooth domains, which yields a new partial data result for Calderoii 
inverse conductivity problem. It is also shown that bisweep data are (up to 
a constant scaling factor) the Schwartz kernel of the relative Neumann-to- 
Dirichlet map. A numerical method for reconstructing the supports of inclu- 
sions from discrete bisweep data is also presented. 



In electrical impedance tomography (EIT) [36] [2] , the conductivity equation 



where D is a bounded domain in or R^, is studied. A fundamental question, the 

so-called Calderon inverse conductivity problem, is whether a measurable conduc- 
tivity a can be determined from the Neumann-to-Dirichlet map (current-to- voltage 
map) Act : / M- u\gD, or equivalently, its inverse, the Dirichlet-to-Neumann map 
A"-'^. In the plane, D C K^, a positive answer was given by Astala and Paiviirinta 
[1] under the assumption that D has a connected complement and a is essentially 
bounded away from zero and infinity (see equation 1.3a below). This followed a 
number of uniqueness proofs in case of more regular conductivities, e.g., [33] [3]. 
The problem has also been extensively studied in dimensions n > 3, see [36] and 
the references therein. 

The conductivity equation is typically studied with boundary currents / € 
H~^/'^{dD), but from a certain point of view it is also natural to consider more sin- 
gular distributional values, in particular, various point distributions on the bound- 
ary of the object. The approach is inspired by the concept of point electrodes [17], 
which can be used as an alternative approximate physics model in place of the so- 
called continuum forward model, the basis of many theoretical results and numerical 
methods in EIT. Another commonly used framework is the com,plete electrode m,odel 
(CEM), which is highly realistic but its theoretical properties remain less studied 



1. Introduction 



(1.1) 



V • (o-Vu) =0 inD 



du 
dv 



f ondD, 



[6][2]. 



We study measurements of the form 



(1.2) 



'w{x,y,p,q) = {Sa:-Sy,{Aa-Ai){6p-dg)), 
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where Ai is the Neumann-to-Dirichlet map for Laplace equation in D. In domains 
with a smooth boundary dD, this type of distributional measurements can be tan- 
gled using the techniques from [29] (cf. [18]). This paper presents an alternative 
approach, based on [27], that is also applicable in piecewise smooth domains. It 
appears that the point measurements w are well-defined provided that 



Thus, in addition to being bounded, a is assumed to be homogeneous near the 
boundary dD. 

In dimension n = 2, there are recent results considering Calderon problem with 
partial data; is it possible to recover a from partial knowledge of A^. (or A~^)? 
It was shown by Imanuvilov, Uhlmann, and Yamamoto [25] that if a is, a priori, 
known to have smoothness of class a > (this assumption is relaxed in [23]), 
then it can be recovered from the knowledge of {{v, {A~^v)\e) ■ suppt; C E} for 
any (relatively) open subset E oia. piecewise smooth boundary dD. A similar result 
for the Neumann-to-Dirichlet map A^ and conductivities a G W^'P{D), p > 2 in 
smooth domains D is formulated in [24]. 

In this paper, we show that if the global smoothness assumptions in [25] [23] [24] 
are exchanged for another, arguably unwanted, restriction (1.3b), a measurable 
conductivity a is uniquely determined by various types of partial data. In particular, 
a result in [21] is generalized by showing that a is recovered from the knowledge 
of ■w{x,y,p,q), x,y G p,q S S for arbitrary countably infinite r,S C dD. In 
addition, a can be recovered from measiu'cments of the type {Ao-f)\w, supp / C V, 
where V,W ^ are arbitrary, possibly disjoint, relatively open subsets of V, W. 

Many numerical reconstruction methods in EIT are easiest to apply if D = B is 
the unit disk. Together with numerical conformal mappings, the concept of hisiueep 
data [21], <;^(x,y) := w{x,y,x,y), offers means of using any such reconstruction 
algorithm in other piecewise smooth, simply connected plane domains. In Section 4, 
we formulate a method for applying unit-disk-bascd reconstruction algorithms in 
other domains, and demonstrate it with the factorization method [4]. 

The main results of this article are outlined in the next section and proven in 
Section 3. For completeness, auxiliary theorems about conformal maps are given in 
an appendix. 



Definition 2.1 (Piecewise smooth domain). Let D (Z he a bounded, simply 
connected Lipschitz domain with a simple boundary dD consisting of a finite number 
m of ^^'"^ smooth arcs (for some < a < 1). We call such D a piecewise "if^'" 
smooth plane domain. 

Throughout this paper, the symbol D is used to denote an arbitrary piecewise 

^i,a smooth domain, which is also identified with the corresponding subset of C 
when necessary. Notice that, as a Lipschitz domain, D cannot have any cusps, 
and the surface (Lebesgue) measure s is well-defined. The symbol v is used for 
the outward unit normal vector of dD, which is defined s-almost everywhere (cf., 
e.g., [14, §1.5]). Sobolev and Lebesgue spaces in D (resp., on dD) are defined with 
respect to the Lebesgue measure in (resp., s). 



(1.3a) 
(1.3b) 



C > a > c> a.e. in D, 
supp(l - cr) CC D. 



2. Setting and main results 
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Let a e L°°{D) satisfy (1.3), which is also assumed in what follows. For any 
/ e H^'^''^{dD), let uf e iJi(£>)/K denote the unique solution to the weak form 
of (1.1) [1][2]: 

(2.1) I aVuf ■Vvdx = {f,jv) ior all v G H'^ {D) /R, 

where 7 : H^{D)/'R H^/^{dD)/'R is the continuous Dirichlet trace operator (sec 
[34, §2.5.4]). Neumann-to-Dirichlet map is the bounded and self-dual operator A^- : 
H^'^^'^idD) H'^/'^{dD)/R defined by / 7^^. We also study the background 
Neumann-to-Dirichlet map, Ai : H^'^''^{dD) H'^/'^{dD)/R, f H> ju{, where 
u{ e H^{D)/M. is the unique solution to (2.1) with (7=1, i.e., Laplace equation 
with Neumann boundary value /. 

By '^'{dD), wc denote the Banach space of all finite real signed Borel measures, 
supported in dD and equipped with the total variation norm (cf. [27]). It is the 
dual space of 'rf{dD) [10, Chapter IV, Theorem 6.2]. The subspace of measures 
fi e 'rf'{dD) such that n{dD) = is denoted 'tfKdD) and it can be identified with 
the dual of '^(9D)/M, the space of continuous functions on dD defined up to a 
constant, which is complete and can be supplied with the norm (cf. [34, §1.1.7]) 

\[fyidD)/R = inf WgyidD) = inf sup |/(a;) + C\. 

The functions in '^oidD) = {/ e 'rf{dD) : /g^/ds = 0} are dense in the 
weak* topology of '^^{dD), that is, for any ^, € '^l{dD), there exists a sequence 
{lJ.j} C %'(i{dD) such that (cf. Lemma 3.14) 

/ flij ds^ [ fdjjL for all / € 'i^{dD)/R, 

JdD JdD 

which is denoted /ij /x. 

Let fi G 'rfl{dD). As a generalization^ of the conductivity equation for measure- 
valued boundary currents, the weak Neumann problem of finding a function e 
W^''^{D)/R = {{u + C : CgR} : w, |Vu| G Li(i:»)} such that 

(2.2) / aVu^" ■V'fidx= [ ipdji for all e '^^(^D) 

J D JdD 

is studied. The corresponding background problem is 

(2.3) / Vw^-V</?dx= / ipdjjL for all V? e '^°°(S), 

J D JdD 



^Unique up to an addition of a constant, that is, in H^{D)/'R := {{u + C : C G K} : ti £ 
i?i(D)} (cf. [34, §1.1.7][30, §1.1.13]). The solution exists if and only if the Neumann boundary 

value has zero mean, that is, / e Ho ^^^{dD) := {g e H-^/'^(dD) : {g, 1) = 0}, which is the dual 
space of H^i''^{dD)/R. The brackets denote dual evaluation. See, e.g., [14] [34], for the definition 
of the space H'-Z'^dD), where it is denoted W^^'^idD). 

hf fi e H-'''\dD)nKidD) that is, M G K(dD) and \Jga'pdfj.\ < C||(p||^i/2(gB)/R. ^ 
which case fj, can be identified with the continuous extension of 93 i-> /g^ (pdn to ^^^[dD), the 
definitions (2.1) and (2.2) coincide. 
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and, as shown in Lemma 3.6, one representative of the background solution is given 

by 

(2.4) <(x)= / 7V($(x),$(j/))dM2/), 

where A'^ is the Neumann-Green function of the unit disk B: 

) if y ^ 0, 

and $:£)—>■ S is a conformal map. 

The main results of this paper are connected to the following concept. 

Theorem 2.1. There exists a unique, linear, and bounded operator, — Ai : 
%''l{dD) ^^{dD)/R, called the relative Neumann-to-Dirichlet map, such that for 

all fx,rie "^lidD), 

(2.6) / [(A,-Ai)M]d7? = g<.(/x,r?) := / (l-a)Vw''-V«Jdx. 

JdD J D 

The map is also self-dual in the sense that the bilinear form Q„ : "^lidD) X 
^l{dD) — > M is symmetric: Qa{t^,rj) — Qcr{f-,v)- Furthermore, 

(2.7) Q„{n, ri) = lim (77^, (A^ - Ai)/Xj), 

where {/Uj}, {rjj} C 'tf^{dD) are any sequences such that jjLj /U and rjj 77. 

The function (A^- — Ai)/i is given by the Dirichlet trace 'y{u^ — u^), where u^^ — 
< e H^{D)/R (LcmmaS.lO). Thus the restriction of A^-Ai to H^^ {dD)r]'^l{dD) 
coincides with the difference of A^. and Ai as defined on page 3, which justifies 
the notation. In smooth domains D, the definition (2.6) likewise coincides with 
the continuous extension Acr — Ai : H^^{dD) (dD)/'R of the relative map to 

distributional Sobolev spaces with s > f /2 (Lemma 3.3). The continuity of A^ — Ai : 
^l{dD) — )• ^(9Z))/R is proven with the aid of a factorization in Theorem 3.8. 

Definition 2.2. The function : dD x dD ^ R, 

'^aix, y) = Qa{5x - 5y, 5x - 5y), 

where 5x f ^ f{x), is called the bisweep data of a. 

Bisweep data can be seen as a point electrode model of the following two-electrode 
EIT measurement: the voltage required to maintain unit current between electrodes 
at X and y is measured. The same measurement is then performed with a homo- 
geneous reference object of shape D (i.e., a = 1 m D). The difference between 
these two measurements, as a function of the electrode positions, is modeled by the 
bisweep data. See [f 7] and [15] for a rigorous treatment of this argument in smooth 
domains. 

Theorem 2.1 allows one to extend the partial data result [21] for Calderon prob- 
lem to pieccwise smooth plane domains: 

Theorem 2.2. LetT C dD be countably infinite and a satisfy (1.3). The knowledge 
of^a on r X r uniquely determines a. 

This can be generalized as follows 



(2.51 



N(x. 11) 



log \y - 



log 



JL 
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Theorem 2.3. LetT,^ be (possibly disjoint) countably infinite subsets ofdD. The 
knowledge of 

(2.8) w{x,y,p,q) := Qa{5x - Sy,Sp - Sg) 

for all x,y gT, p,q uniquely determines a. 

The above also relates to the next corollary, which can alternatively be proven 
directly using a unique continuation argument. (Cauchy data for Laplace's equation 
in some neighborhood of the boundary is always available on a certain part of dD.) 

Corollary 2.4. Let V,W be (possibly disjoint) non-empty relatively open subsets 

of dD and denote by '^o(y) C '>ao{dD) (resp., '^o{W)) the continuous mean-free 
functions supported in V (resp., W). The knowledge of 

{if,gAf,K9)) ■■ fG'^oiv), g&'rf.iw)}, 

uniquely determines any conductivity a that satisfies (1.3). In other words, if{f, A.„g) 
= (/, A^g) for all f G ^o(^)> 9 € '^o(W^); then a = a almost everywhere in D. 

These results are based on the fact that, in the unit disk B = D, bisweep data 
and its generalization in Theorem 2.3 are jointly analytic functions. Bisweep data 
also comprise (up to a scaling factor) the Schwartz kernel of the relative Neumann- 
to-Dirichlet map in the sense that 

Theorem 2.5. For any /i,r/ e %{dD), 
/ [(A^ - Ai)/i]dr? 

JdD 

and <;„ € '^°'{dD x dD) for some a > 

Many of these relations can be proven by generalizing the corresponding result 
from the unit disk B (or some other smooth domain) to an arbitrary piecewise 
smooth plane domain D using the fact that point current sources /j, = YyjLi o^j^xj 
are moved naturally by conformal mappings (cf. [15]). In particular, bisweep data 
provide a natural method of "transporting" relative Neumann-to-Dirichlet maps 
between different domains. Namely, if a = (7o$~^, where $ is a conformal map 
from D to the unit disk B, then 

(2.9) <;,ix,y)^^^iMx),Hy)) 

for all a;, J/ e dD. Notice that, since dD is a Jordan curve, $ extends to a homeomor- 
phism between D and B. In fact, $ is Holder continuous and its further smoothness 
properties are determined by the smoothness of the arcs of dD and its vertex angles 
(cf. [35, Ch. 3] and Theorem A.2). 

3. Proofs of the results 

It follows from Poincare inequality and the boundedness assumptions on a that 
the left hand side of (2.1) defines a bounded and coercive bilinear form in the Hilbert 
space H^{D)/M.. The right hand side is a continuous functional in {H^{D)/M.y due 
to the trace theorem [34, §2.5.4]. Hence the unique solvability of (2.1) and continuity 
with respect to the data follows readily from Riesz representation theorem, (cf. [8, 
Chapter VII, §1.2.2]) 

In case of (2.2) and (2.3) the above technique fails because the distributions 
defined by the right hand sides are not generally bounded in (iJ^ (£)))' (since 



JJ <;„diJLXri, 

dDxdD 
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%"{dD) ^ H~^^^{dD)). To the best knowledge of the author, the distributional 
theory studied by, e.g., Lions, Magenes, and Necas [29] [34] is not directly applica- 
ble either, due to the higher regularity requirements for the domain D. 

This section shows an alternative approach, based on the work of Krai [27], for 
solving these problems. 

3.1. Weak Neumann problems. The fact that the weak background Neumann 

problem (2.3) has a unique (harmonic) solution ?ij G '«f°°(-D)/IR such that jVu^j € 
L^{D) follows from [27] provided that certain geometric assumptions on D are 
satisfied. It is stated in, e.g., [32] that piecewise "^^'^ smoothness is sufficient. 
Regarding the conductivity equation, there holds 

Lemma 3.1. The weak conductivity equation (2.2) has a unique solution u'^ € 
Wi'i(£))/M, given by = u;^, where G H^{D)/R satisfies 



(3 



.1) [ aVw'' ■\'ipdx= [ {1- a)Vu'^ -Vipdx for all if G H'^ (D) /R. 

J D J D 



Proof. Due to Poincare inequality [34, §1.1.7] and the boundedness assumptions 

in (1.3), the left hand side of (3.1) is bounded and coercive. The right hand side 
defines a continuous functional ^p ^ Jjj{l — a)\/ui-'V(pdx in (i?^(£))/M)', for is 
smooth in supp(l — a). As a result, there exists a unique solution w^. Furthermore, 
'r°°(D)/M is dense in H^{D)/m. (cf. [30, §1.1.6]) so (3.1) is satisfied for all (p G 
H\D)/R if and only if it holds for all (p G 'rf°°(D)/R. It follows that + w*" is 
the unique solution to (2.2). 

Notice that < is also in V'{D)/R by [30, §1.1.11] and therefore <, w'' G 
W^'^{D)/R. □ 

The following two lemmas state the relationship between the weak formulations 
presented above and a distributional Sobolev space formulation based on trace 
theorems, utilized in, e.g., [21]. 

Lemma 3.2. Let D he a smooth domain, fi G 'rfl{dD) arbitrary, and s < 1. 
Then Ui = solves (2.3) if and only if Ui G and 

(3.2) Aui =0 in D, 'fiUi = ^ on dD, 

where 71 : X** — ^ H'^~^/'^{dD) is the continuous extension of the normal derivative 
in, e.g., the Banach space = {u & H^{D)/R : Au G i^(D)} equipped with the 
graph norm (cf. [29, Chapter 2, Remark 7.2]^. 

Proof. Assume that = solves (2.3) and s < —2. Clearly, Aui = in D. 
According to [30, Theorem 1.1.6.2] (and its proof), there exists a sequence {uj} C 
'^°°{D) that converges to (an arbitrary representative of) w.r.t. the norm 

ll^ll* hWmD) + |||Vw|||Li(r,) + ||A-(;||i2(£,). 

Due to Sobolev's embedding theorem. 



= sup 



L 



vip dx 

D 



< 



sup II II LI (D) II ¥'11 ^(5) < C\\v\\l^d) 



=1 



and consequently, uj — )■ wi G X ^, — )• 71U1 in H '^/^{dD). For any ip G 
^°°{dD), let (p be an arbitrary extension of ip to '^°°{D). Then, by Green's first 
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identity, 

(7iUi,<^)= lim(^,(^)= lim / ^(y?ds= lim / (Vu,- • V(p + (pAu,) da; 
= j Vui • V<^ dx = ip d/x, 

Jd JdD 

which means that 71U1 = /x. Thus Ui solves (3.2). 

Since the solution Ui to (3.2) with s < — 2 is unique, it follows that it can always 
be identified with a function (equivalence class) ui that satisfies (2.3). If ui solves 

(3.2) for —2 < s < 1, this also remains true for s < —2. Furthermore, a unique 
solution to (3.2) exists for any s < 1 (because '^'{dD) C H-'^/^-''{dD) for any 
e > 0), which proves the general claim. □ 

In smooth domains D, it is possible to extend the operator Ao-— Ai : H^, (dD) — > 
H^/^{dD)/R (as defined on page 3) to a continuous mapping between the Sobolcv 
spaces H-^idD) and H-'{dD)/R for any s e M (cf. [18]). The next lemma shows 
that the result coincides with the definition (2.6). 

Lemma 3.3. Let D be a smooth domain and s > 1/2. Then 

(3.3) {t], {A^ - Ai )/i) = (/i, T]) = Q„ (?7, n) 

where ii,r] G '^l{dD) are arbitrary and A„ - Ai : H-^{dD) H^{dD)/R. 

Proof. Fix an arbitrary s > 1/2 and let U CC D be such that supp(l — a) CC U. 

For any <j> E H^^{dD), lot ?if be the unique solution to (3.2) with boimdary value 
(/), which is clearly equivalent to definitions (2.1) and (2.3) for smooth 0. Moreover, 
for any smooth functions ^, </3 € '^^{dD), the "Green's formulas" 



(3.4) 



{if, {A^ - Ai)<P) = / (1 - cr)Vut ■ Vuf dx= I (1 - c7)Vu'^ • V< da; 
J D Jd 

= [ (l-£7)V(4 + w^)-V<da; 
Ju 



where w'^ solves (3.1), follow readily from the definitions of A„, Ai on page 3 and 
Lemma 3.1. 

Due to standard elliptic regularity theory (cf., e.g., [18, Appendix]), the operator 
(/> H> w^|(7 : H-^idD) H^{U)/R is continuous. The mapping uf\u H> w'^ : 
H'^{U)/R H'^{D)/R is also bounded (see the proof of Theorem 3.8 for details). 
Consequently, (3.4) is well-defined for any 4>, (fi G H~'^{dD) and it yields the unique 
continuous extension A^ — Ai : H~^{dD) — >■ H^(dD)/R. It follows from Lemma 3.2 
that the extension satisfies (3.3) for any /i, € ^l{dD) c H~^{dD). □ 

In the next section, the formulas (2.2), (2.3) need to be considered with less 
smooth test functions ip. The following lemma states that they remain valid as long 
as is Lipschitz. 

Lemma 3.4. If u = solves (2.2) then (2.2) in fact holds for all G "^"'^(L)). 

Proof. Since (f is Lipschitz continuous, there exists an extension cp : M'' — > M such 
that (p\jy — p and || |Vi^| [jioo^i^d) < C (cf. [11, §5.8.2] and [31]). Let ipc = * ip G 
where G 'ifQ°{B^{0)) is a mollifier and * denotes convolution [11, §C.4]. 
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Denote Djj = {x E D : dist{x, dD) > rj} (where i] > e). Then 

/ ctVm • Vife dx — aVu ■ V<^ dx = I aVu • (Vi^^ — \7ip) dx 

J D J D J D 

< IkllL-(D) ( 2C /" \Vu\dx+ [ \Vu\\V^-V^,\dx\, 
\ Jd\d^ Jd„ J 

which can be made arbitrarily small by first choosing a suitable 77 > and then 
e < T], because — ^ G ^^i^ri) and Vu € L^(D^). On the other hand, 

/ aVu ■ Vipe dx = iff: d/U ^~^°> / i^dji 
Jd JdD JdD 

due to the uniform continuity of ip. □ 

3.2. Conformally mapped Neumann problems. A key ingredient in the par- 
tial data results in Theorems 2.2 and 2.3 is the ability to transform to an equivalent 
problem in the unit disk. To this end, the next theorem describes how the solution 
to the weak conductivity equation transforms under conformal mappings between 
piecewise smooth domains. 

Theorem 3.5. Let D,E be piecewise smooth plane domains and ^ : D ^ E 
a conformal map. If iif^ solves (2.2), then m = w,^ o solves 

(3.5) / aS7u-S7>fidx= [ <pdfi for all ip e '^'^ (E) , 

Je J as 

where 

JdD 

Proof. Let us first study the case when $ e ^^(D) (i.e., $ can be extended to 
such a function). For any (p e '^°°{E), let (p = (p o ^. Then tp € ^^(Z?) and, by 
Lemma 3.4, 

/ (pdjl= I <^d/i= / ctVu'' • V<^ dx = / aVu-S/ipdx, 
Joe JdD Jd Je 

where the last step follows from Lemma A. 3. 

Conversely, if u solves (3.5), then it must hold that u — o (f>^^, because the 
solution to (3.5) is known to be unique (up to addition of a constant function). 
Thus u o $ solves (2.2). 

Due to Theorem A. 2, there exists, for some a' > 0, piecewise "if^'" smooth 
domains E' and conformal maps $f , $f , $f , $f € such that 

D^D'^B^E'^E 

and $ = ^2 ° (^f ° ^? ° (^f)"^- The claimed mapping property is known to 
hold between any pair of consecutive domains in the above chain since either the 
relevant conformal map or its inverse is smooth enough. This proves the general 
claim. □ 

It follows immediately from the above theorem that = o $ and w^^ = 
u'^ — Ui = o $ transform similarly. Also observe that point current sources are 
not "deformed" by the transformation; if = CjSx^, then fi = J2j Cj^*(xj)- 
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Lemma 3.6. Let ^ : D ^ B be a conformal map. The representation formula 
(2.4) holds and, furthermore, 

(3.6) V<(x) = / V. (Ni^x), i.(y))) d/x(t/) 

JdD 

for any x € D. 

Proof. Let U CC B be arbitrary. The explicit representation (2.5) shows that 
^^Pxeu,y(^dB \^{x,y)\,sup^^uy^gB \'^xN{x,y)\ < oo and it is known that the op- 
erator / r4\u ■■ H-'{dB) H^{U)/R, s eR, where u{ = u{ solves (3.2) in 
D = B,is bounded [18, Appendix]. It follows by a straightforward density argument 
that 

(3.7) {N{x,-),f), xGU, 

defines a representative of u{\u and {'VxN{x, •), /), x G U, defines Vu{\u- Since U 
was arbitrary, the above expressions remain valid formulas for (a representative of) 
u{ and Vu^ in the whole disk B. 

The general claim for an arbitrary piecewise "^^'^ smooth D now follows from 
Lemma 3.2 and Theorem 3.5, since 

< = <o$=/' iV($(.),y)d/i(y)= / 7V($(.),$(y))d^(2/) 

Job JdD 

modulo constant functions. □ 

Remark 3.7. The representative given by (3.7) (for / e ^{dB)) has zero mean on 
the unit circle Jg^ u{ ds = 0. Correspondingly, the representative given by (2.4) 
(for any g € '^{dD)) satisfies 

/ ulo $-1 ds = 0, 

JdB 

and its mean does not generally vanish on dD. Thus N{<t{-), $(•)) is not the usual 
Neumann-Green function of D (cf., e.g., [4]), but induces a different normalization 
criterion ("ground level") for u^. 

The rest of this section focuses on proving the claimed properties of — Ai in 
Theorem 2.1. 

Theorem 3.8 (Factorization). Let U CC D be an open set such that supp{l — a) C 
U. Then for any iJ,,r] e '^l{dD), 

(3.8a) Q„{n, r])= f {SGfi) ■ {Gr,) dx, 

Ju 

(3.8b) = / {l-a)\/{A7])-S/{{I + P)Ai^}dx, 

Jd 

where the operators 

G:%i(dD)^{L\U)f, M^V<|c/, 

P : H^{U)/R^ H^{D)IR, u'i\u ^ w^" , 

and S : (L^(C/))^ — > (L^([/))^ are linear and bounded. In addition, G is continuous 
between the weak* topology of^KdD) and the strong topology of {L^{U))^. 
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Proof. The continuity of the (well-defined) operators G and A follow from Lem- 
ma 3.6: 



\iAn){x)\ 



[ iV($(x),$(y))d/x(y) 

JdD 



< sup \N{x,y)\-\\iJ,yngD), 

£e*(!7), yedB 



which is finite, as seen from the explicit formula (2.5). Similarly, |(G/i)(a;)| = 
\Vx{A.iJ,){x)\ < C||/u||<^^(a£)). Furthermore, if /z, then {Giij){x) converges 

pointwise for x G U by (3.6) and 

\{G^,j){x)\< sup |V.iV($(ar),$(2/))|-||M,|k'(aD)<C, 

since since a weak* convergent sequence is strongly bounded. Therefore \ {GiJ,j){x) — 
(G/i)(.T)p is uniformly bounded for x G J7 and it follows from the dominated con- 
vergence theorem that G^j — >■ Gji in (L^([/))^. 
Now consider the variational problem 



aVw ■S/vdx = f{v) for all v e (D) /R. 

ID 

As in the proof of Lemma 3.1, for any / e iH^{D)/Ry, there exists a unique 
solution w e H^{D)/R and the mapping T : {H^{D)/RY H^{D)/R,f ^ w is 
continuous. 

Set F : H^{U)/R {H^{D)/Ry, g ^ {v ^ - g -Vv ^) and P = TF: 
H\U)/R -J> H\D)/R, both of which are bounded. Then 

= T[v^ /[/(I - cr)V< • Vvdx] = TFu^lu = PAji. 
It follows from definition (2.6) and Lemma 3.1 that 

Qa{^^, rt)= /" (1 - <^)V« + w^) ■ Vul dx 
Ju 

(1 - ct)V((7 + P)Aij) ■ V{Ari) dx. 



Similarly, set F : {L^{U)f {H\D)/RY, g ^ {v ^ J^{l-a)g-yvdx), whence 

w^" = TFVul = TFG 
and the factorization (3.8a) holds for 5 = (1 - cr){I + VTF). □ 
Lemma 3.9. The bilinear form Q„ is symmetric. 

Proof. Let $ : D — )■ i? be a conformal map. Due to Theorem 3.5 and integration 
by substitution (cf. equation A.l), 

Qa{t^,ri)= [ {l-a)Vu>'-Vu1dx= [ (1 ct)Vu'^ • Vwf dx 

J D J B 

where a = a o and p,, fj are the weak Neumann boundary values of o 
and ul o respectively. By Lemma 3.3, this equals 

( (1 - a)V< • Vw^ dx = [ (l-a)V<- VM''dx = Q<,(r/,At). □ 

J B J D 

So far, we have shown that is well-defined and bounded. However, it remains 
to prove that, for each /x S 'r^l{dD), there exists a unique (A^^ — Ai)// G 
such that (2.6) is satisfied. 
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Lemma 3.10 (Continuity of (A^ — Ai)/^). Let G Wl{dD) and be as in (3.1). 
Then {A^-Ai)fi = e '^{dD)/R. 

Proof. Let jl, and ct be as in Theorem 3.5. Since w'^ (as in equation 3.1 with 

D = B and a = a), is harmonic near the boundary dB and satisfies a homogeneous 
Neumann condition, w'^lj^^fj G 'if°°(B \ U)/M. for some U CC B. By Theorem 3.5, 
— o ^ where $ € '^{D). Therefore the trace coincides with the pointwise 
Umit of a uniformly continuous function (equivalence class) w^\jj\^^^r^ on the 
boundary dD and 7^^ G '^(dD) /R. 

For an arbitrary / G "^0(31)) C H^^''^{dD), 

I -fwi^fds = ifnw^) = [ o-Vw'' • Vu^ dx 

JdD J D 

= f {l-a)Vu1-VuUx = QM,l^)=QAf^J)- 
Jd 

It now follows from Theorem 3.8 and the weak* density of ^o(9£') in '^^{dD) that 
for any r] G %{dD), 

j jw^dr]= lim / ^w^rijds= lim Qcr{lJ',ilj) = QaifJ-iV)- 

Evaluating /^^[(Act — Ai)/i] d{Sx—Sy) for all x,y G dD also shows that (A^ — Ai)/x : 
dD — >• M is unique up to a constant. □ 



The continuity of A^- - Ai : "^UdD) '^{dD)/R, as well as (2.7), follow from 
Theorem 3.8, Lemma 3.9, and the boundedness of weak* convergent sequences. 



3.3. Bisweep data. It is known by Lemma 3.3 that, in the unit disk B, Defini- 
tion 2.2 coincides with the definition of bisweep data in [21]. The identity (2.9) 
follows directly from Theorem 3.5 and thus we are ready to state: 



Proof of Theorem 2.2. Let $ be a conformal map from D to the unit disk B and 
f = $(r). By [21, Corollary 2.3 & Remark 2.4], the knowledge of ^{x,y) for all 
{x,y) G TxT determines the unit disk conductivity a and a is recovered as cro$. □ 

The cornerstones of these results are the solvability of Calderon problem with 
full data [1] and the analyticity of bisweep data, which is also proven in [21]. For 
completeness, wc will give a (slightly extended) proof of this fact, which will also be 
utilized in the proofs of Theorems 2.3 and 2.5. Here the "interior domain" factor- 
ization (3.8a) is utilized instead of the transmission-problem-based factorizations 
of the Neumann-to-Dirichlet map in [15], [22], [21] and [18]. 

Theorem 3.11. Let D = B be the unit disk and 6{6) := ^(cose.sine)- '^he function 

is (jointly) analytic in K^. Furthermore, it extends to a complex analytic function 
in C C*, where V = {zgC: \^z\ < X) for some I > 0. 
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Proof. Consider the function 



(3.9) 



1 

2^ 




where = Xi + ix2- The ultimate step on the first fine follows from a straightfor- 
ward complexification of (3.6) in the case $ = Id. 

Let < i? < 1 be such that supp(l - a) CC U := Br{0). Set I = -^logi? 
and let V be as above. It is clear from the explicit representation (3.9) that each 
component of g{x, •, •) extends to a holomorphic function of 6,(p gV when x G U. 
In particular, there exists an m > such that — e*^| > m and — e~*^| > m for 
all 9 gV, X £ U. Hence the real and imaginary parts of g{-,9, </>) are in (L^([/))^. 

The real-valued operator S in (3.8a) can be naturally complexified as S{u+iv) = 
Su + iSv. Thus the extension of 



(3.10) 



o^i^i, O2, 03, 64) 



Ju 



e2)-iSyg{y,es,e4)){x)dx 



to is well-defined. 

It it apparent from the explicit representation (3.9) that for any 6,(1) G V and 
7] gC\ {0} such that 9 + r] G V, there exists a bound M > so that the inequality 



g{x,9 + r],(f))-g{x,9,(f>) _ dg 

rj 89^ ' '^^ 



< M\r}\ 



holds for all X G t/. As a result, for any bounded linear mapping L 
{L^{U))^ (in particular, Id and S), 



{L\U)f 



{Lg{;9 + n,<l>))-{Lg{;9,<j>)) 



V 

g{-,9 + v,(l})- 



89 



< \\L\\M\rj\ 



LHU) 

The same holds true for the other complex variable (j) of g. This means that the 
operator {9, (p) ^ Lg{-,9, (p) : V"^ ^ {L'^{B)Y is strongly holomorphic in both vari- 
ables. (The complex partial derivative w.r.t. 6 is (^,^) Z/(||(-, 0, ^)).) Therefore 
the expression (3.10), which is a continuous bilinear form in (L^(C/))^, is separately 
holomorphic in V with respect to each variable while the others assume arbitrary 
fixed values in V . The claim of joint holomorphy in follows from Hartog's The- 
orem [20, Theorem 2.2.8], and the analjd:icity in follows via restriction. □ 

By similar arguments (cf. [21]), one can prove the following 
Corollary 3.12. The mapping w : {dB^ K^, 

^(2^1, ^2,^^3,2^4) = Qai^z-, - ^;J2>^^3 - ^Zi) 

extends to a complex analytic function in some neighbourhood of [dB^ C C^. 
Together with (2.7), the next lemma proves Theorem 2.5. 
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Lemma 3.13. For any f,g E '^^^{dD), it holds that 

Q<T{f,g) = -\ jj f{x)g{y)<;{x,y)ds{x)ds{y). 

ODxdD 

In addition, 

N 

(3.11) <;^{x,y)= lim ^^Xik {(t>k{x)(t>i{x) - 2(t>k{x)(t>i{y) + 4>k{y)(t>i{y)) 

i,k=l 

uniformly in dD x dD, where is a certain orthonormal basis for L'^{dD) 

and Xij = Qa{<j)i,(j)j), i,j = 1,2,3, .. . 

Proof. Let cq € ^"([0, L)) be a natural parametrization of dD, i.e., f^^f^f^^ ds = 
& — a for all < a < 6 < i. Denote by c the periodic continuation of Cq to ^"(M) 
and let 

Wj{x)}^^o = { ^, ^1 COS ^f,^sm^^,^ cos . . .| 

be the standard trigonometric Fourier basis of L^{{0,L)). 

Let $ : — ^ B be a conformal map and a = a o Corollary 3.12 shows 

that the function w{x, y,p, q) — Qai^x — 5y,5p — 5q) is smooth (jointly analytic) on 
dB X dB X dB X dB. Thus 

^{x' ,y' ,p' := Qa{S^{c{x')) - ^3>(c(y'))'^*(c(p')) ~ ^*(c(g'))) 

is Holder contimious in M.^ (because $ G "^"(D)) and therefore the following four- 
dimensional Fourier series of ui, where ifiijki{x' , y',p', q') = (pi{x')(fj{y')(pk{p')'fii{q'), 
converges uniformly [13]: 

N 



w = lim (pijki / (fijkii^ dt^ 

N—^oo ^ — ' / 



Now set (f>j = Lpj o . By Theorem 3.5, Lemma 3.10, and definition (2.8), 

w{x,y,p,q) = uj{cQ^{x),c^^{y),c^^{p),c^^{q)) = w'P^'^ix) - w'^-'^iy), 

where w^''^ e '^<^{dD) denotes the continuous representative of the Dirichlet trace 
jui]Sp-Sg -^j-tj^ 2ero mean and w^^~^'' is as in (3.1). 

Similarly, let Wi denote the representative of ^w"^^ with zero mean. By Theo- 
rem 2.1 and Lemma 3.10, /^^ wP''^4'i ds = Wi{p)—Wi{q) and \ij = Xji = Jgj^ Wicpj ds 
for all i,j = 0,1,2,... if one sets Aq, = Ajo = 0. Define (l)ij,4'ijk,4>ijki like (pijki 
above. Then 



/ ipijkiojdt^= / (j>ijkiw ds'^ 

= [ ^jki [ {wP''^{x)-wP'\y))Mx)ds{x)ds' 
J{aD)3 JdD 

= [ ct>jki{l w^'Vids- W(y)iO 

J(dD)^ \JdD J 

= / (t>jki ({wi{p) - M<l)) - SiowP''^{y)Li) ds^ 
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= / ^fci {Sjo{wi(p) - Wi{q)) - Sioiwjip) - Wj{q))) ds^ 
= 4>i [SjoiXik - hoL^Wi{q)) - SioiXjk - hoL^Wjiq))) ds 

JdD ^ ' 

= {SjoiSioXik — (5fcoAj|) — 5io{5ioXjk — Sko^ji)) L, 

where dij is the Kronecker delta. 
Inserting the above expression into 

N 

'iaix, y) = w{x,y,x,y) = \im V <pi{x)(j)j{y)<l)k{x)(pi(y) / fijkiujdt'^ 

iv— >oo ' J 

i,j,k,l=0 (o_^)4 

and simplifying yields (3.11) and the limit is uniform in dD x dD. 
As a result, 

N 



#^Ay^z)f{y)g{z)ds{y)ds{z) = -2 hm V A^j / ^Jds / ^jgds 
J9D JdD 

= -2 Jim^ J^^ (E^Ii loD "t^if ((^-^ " ) (Ef=i '/'i /an ^^'jfi' ds) ) ds 
= -2[ /((A,-Ai)5)ds, 

where the last equality follows from the boundedness of A^- Ai : L^{dD) — > 
L^{dD)/R and the fact that E»=i 0i /od "^^i/ ds ^ / in L^{dD) for any / e 
Ll{dD). □ 

The analyticity of the "four-electrode function" Co can be combined with the 
conformal mapping property of point current patterns, which yields the new partial 
data result: 

Proof of Theorem 2.3. Let $ : ^ i? be a conformal map. For an arbitrary branch 
of the complex logarithm, set f = -ilog$(r) C M, S = -■jlog<I>(S) C R. By 
Theorem 3.5, the measurements (2.8) determine the quantity Cj{6, Oq, (p, (po), defined 
in Theorem 3.11, for all 6, Oq G F, ^, (pQ € S. 

Let G r, '/'o G S be arbitrary and ^i, (pi be some accumulation points of T 
and S, respectively. As in [21], for any integers j, fc > 1, it holds that 

^^^(ei,^o,<^i,<^o) = (4f ,(A^- AO^g^ 

where : / ^ ^f{d) and Aj- Ai : H-^{dB) H''{dB)/R is as in Lemma 3.3 
for some sufficiently large s. Furthermore, the values of the partial derivatives on 

the left hand side arc revealed by the measurements. 

Due to the analyticity of {A^ — Ai)(5^^'' e '^^°°(9i3)/M, this also determines 
(A^ - Ai)^^*;^) for all j,k > 1. By [21, Theorem 2.2], Theorem 2.5 and (2.9), 

this yields A^ - Ai : H^^''^{dD) H'^/'^{dD)/R. Since Ai is known a priori, one 
recovers a G L°°{D) from A^. by [1]. □ 
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Lemma 3.14 (Approximation of measures by continuous functions). Let ji G 
'^^{dD) he arbitrary and V any relatively open subset ofdD such that supp CCV. 
Then there exists a sequence {fij} C '^o{V) such that fi. 

Proof. Using a paramctrization of dD and convolution with mollificrs in K (cf., 
e.g., [30, §1.1.5][34, §2.1.3]), one can construct a sequence {jlj} C 'rf{dD) such that 
jlj fi. Let av € '^{V) be such that supp(ay — 1) fl supp/x = and set Cj = 
— fgjj ayjij ds/ fgjj av ds. Then /Xj = [jlj + Cj)av has the claimed properties. □ 

The measurements (2.8) can be interpreted as follows: Let x,y,p,q G dD be 
arbitrary. There exists sequences {,fj},{gj} C 'lo^idD) such that fj 5x — 5y, 
9j ~ ^q- to Theorem 3.5 and (2.7), 

u!{x,y,p,q)^ lim (/,-, (A„ - Ai)^^), 

j->-oo 

where 

if J, {K - ^i)9j) = {fj,Kgj) - (/j, Aiffj), 
and {fj,Aigj) does not depend on a but only the shape of D and is thus known 
a priori. Notice also that, by Lemma 3.14, the sequences may be chosen so that 
{fj} C %{V), {gj} C '^o(H^), which yields Corollary 2.4. 

4. Reconstruction of a from bisweep data 

4.1. Reconstruction method. Given a conformal map ^ : D ^ B one can 

compute a reconstruction of the unit disk conductivity a = cr o from the 
bisweep data and map it back to D using Due to (2.9), c,^ can he computed 
from Co-- In addition. Theorem 2.5 provides means for converting between bisweep 
data and the relative Neumann-to-Dirichlet map A^ — Ai : Ll{D) — > L^(I?)/R. 
However, there are some subtleties concerning discretization and relationship to 
electrode measurements in this approach. 

We consider the setting where ^ : D ^ B is a. conformal map, Xi,. . . ,Xn € dB 
equispaced points, and the bisweep data = are given for all i,j 

I :■; . where xj := ^~^{xj). Prom a practical point of view, this approximates a 

measurement where one can select the positions xi, . . . ,Xn of n electrodes on the 
boundary dD of the object according to a pre-defined pattern, and then conduct 
the (relative) EIT measurements (cf. [17]). However, unlike in smooth domains, 
in piecewise smooth D, this relationship between relative point electrode measure- 
ments and the rcialistic complete electrode model [6] has not been rigorously studied. 

By (2.9), the discrete bisweep data Qj correspond to the unit disk data ?g- sampled 
on the regular grid x {xi, . . . ,Xn} and, due to Theorem 2.5, their 

discrete two-dimensional Fourier coefficients 

^ij ■= X! —4>i{xi)(t>j{xk)<ikU 1 < «,i < n - 1 

kl=l 

with respect to (for example) the basis 

{4>i{{cos6,sm6))}'^I^ = |^,yicos0,yisin0,yicos(20),... 
can be used as an approximation of 



(4.1) 



\j = Qa{<l>i,<i>j) = {<l>j, (Aa - Ai)(^j). 
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Now one can compute a reconstruction of the conductivity a in arbitrary points 
Vit ■ ■ iVm € D by reconstructing a in $(y i ),..., from the approximate 

matrix representation [Xij\i.j=i^....n-i of the relative Neumann-to-Dirichlet map for 
the unit disk conductivity a. In the foUowing numerical examples, the unit disk 
reconstruction is done using the factorization method (see Section 4.2 below), which 
aims to recover the support of the inclusions 1 — a. 

Remark 4.1. If the discrete measurement on dD does not correspond to bisweep 
data but is (interpreted as), for example, Ajj, i,j = 1, ... ,n — 1 (cf. eq. 4.1), then 
one could compute an approximation of Cij by truncating the series (3.11). 

If the measurement can be interpreted as a general measurement with n point 
electrodes, that is, 

(n n \ 

k=l 1=1 ) 

where ai, . . . , a„_i € M" and /3i, . . . , /3„_i € M" are linearly independent sets of 
mean-free vectors, then one can compute Qj from mjj by solving a system of linear 
equations. 

4.2. Factorization method. In this section, a numerical method for locating the 
inhomogeneities in ct e L^{B) (similar to that in, e.g., [19] and [28]) is presented. 
A theoretical result behind the method — originally from Kirsch [26] — is stated for 
completeness, but all further analysis is omitted. The applicability is demonstrated 
by the numerical examples in Section 4.3. For proofs and other properties of the 
factorization method, refer to, e.g., [4, 5, 16, 12, 28, 19]. 

Assume fl CC B is such that B\Q is connected and a G L°°{B) satisfies 

(4.2) 0<c<a<C <1 a.e. in = 1 in B \ H. 

Let {ffe} C Ll{dB) be orthonormal eigenfunctions and {Afc} C M the corresponding 
eigenvalues of the compact operator — Ai : Ll{dB) L'^{dB)/R: 

(As - Ai)Dfc = AfcUfe +R, fc = 1,2,3,... 

Denote c/^,d :dB^R, [4] 

/X , „ ,w X 1 (z — x) ■ d 
where z e B, d € dB cM."^ and N is given by (2.5). 

Theorem 4.2 (Factorization method). Let d € dB be arbitrary. For any z e 

B \ d^l, z £ if and only if [12] 



(4-3) f'^^'^--=WA-^ ^ 

\\gz.d\\L2{aB) k=i 



\JoB 9z,dVk ds|^ 
\\k\ 



< oo. 



This is applied as the following algorithm for reconstructing the inclusion ^l from 
the (approximate) measurements (4.1). First, choose a reconstruction order M < n 
and compute the singular value decomposition 



Lvk = (TkUk, L^Uk = (TkVk, k = l,...,M 
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of the M X M matrix L = [Xij]fj^i that depicts the relative Neumann-to-Dirichlet 
map Ag. — Ai in a finite trigonometric basis. The following truncation 

f , . _ Efcli \g^vk\^/Wk\ 

where g ~ [gi, . . . , ^m] is a vector of Fourier coefficients of gz,d in the same basis, is 
used to "approximate" (4.3). The reconstruction is given hy = {z G B : Ind(2;) < 
Coo} for some cut-off value Coo > of the function 

(4.4) Ind(z) = ^4(z)/iVd 

fc=i 

where some (odd) number Nj, of equispaced points dk, k = 1, . . . , are used as 
the dipolc directions to reduce artefacts (cf. [28]). In practice, Ind is sampled on 
some finite grid (and with different values of Ceo) to produce an image of 51. 

4.3. Numerical examples. In the following examples, the domain D is a polygon, 
which contains an inhomogeneity defined by 



a{x) 



if a; e fii 
otherwise, 



where < < 1 and the simply connected inclusions flj CC D, i = 1, . . . , M are 
strictly separated, whence (1.3) and (4.2) are satisfied. 

Approximate bisweep data are computed numerically from (2.4) and (3.8b), 
where the operators A and P have been discretized using the finite element method 
(FEM). In more detail, the simulation is carried out as follows: 

(1) Construct a triangular finite element mesh (e.g.. Figure la) with nodes 
Zi,. . . ,zn € D, and a corresponding piecewise linear basis C ^^{D) 
such that each basis function (pj is supported on the triangles that are 
adjacent to mesh node zj. Define 

In = {j : supp{ipj) n O 7^ 0} C {1, . . . , A^}, 

and U = Uj£/f5 supp(<^j). Also compute the stiffness matrices 

{K^)ij= a\/(pt{x) ■\/ipj{x)dx, i,j = l,...,N 

J D 

and G M^><^. 

(2) For all x,y € {x\, . . . , a;„} C dD, compute 

u^.^y = A{5^ - 5y){z^) = Ni^zj), ^x)) - Ni^zj), $(y)) 
for each j e In and set Uj'^ = for all j ^ In- This defines a finite 
element approximation u^'^ = Ei^i '"Vj of ^{^x — Sy). Denote Ux,y = 

(3) A finite element approximation w^'^ = EjLi "^j'^^j of Pu^'^ can be solved 
from the linear system 

and w ^{x,y), x = Xi, y = Xj can be computed from (3.8b) as 
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In the examples below, these computations were done in MATLAB and the confor- 
mal map $ was computed using the Schwartz-Christoffel toolbox [9]. 

Remark 4.3. It is also possible to simulate the data by first computing a = cro$~i, 
then computing as described in, for example, [15] and finally mapping the result 
back to D in order to obtain c^. This approach, which is also the theoretical basis of 
the reconstruction method, is not taken in order to avoid an obvious inverse crime 
[7]. In the chosen alternative procedure, the numerical conformal map $ that is 
used for reconstruction is employed only in the simulation of Ui\u- 

In the following numerical examples, the domain D is a non-convex polygon 

depicted in Figures 1 and 2. There arc two inclusions, a disk and a rectangle, which 
have the conductivity ki = K2 = |. A conformal map ^ : D ^ B is constructed 
using the Schwartz-Christoffel toolbox and n point electrodes are positioned as 
Xj = ^^^{xj), where Xj are equispaced on dB. We use Nj, = 15 dipole directions 
for (4.4) in both examples. 

Example 4.4. In the first simulation, a high number n = 128 of point electrodes is 
used. The positions of the inclusions and the electrodes are shown in Figure Id. The 

unit disk phantom a = ao^~^ g L°°{B) and xj are depicted in Figure Ic. Discrete 
bisweep data i^ij = ^a{xi, Xj), i,j = 1, . . . ,n (shown in Figure lb), are computed as 
described above. 

No artificial noise is added to the simulation and the reconstruction is conducted 
with a high order of M = 64 singular values. The result is also shown in Figure Id, 
where the different shades of gray correspond to different cut-off levels Coo (loga- 
rithmic spacing). This suggests that in an ideal setting, the method works as desired 
and the supports of the inclusions are recovered well. 

Example 4.5. In the second example, the number of electrodes is reduced to n = 
16 and artificial errors are generated as follows: The vertices of the domain are 
perturbed slightly, which yields a new domain D' . Perturbed electrode positions 
Xj e dD' are computed by selecting the closest point to Xj on partial D' and adding 
a small perturbation. Bisweep data ^-j = <;„'{x'^,Xj) are computed and extra noise 
(normally distributed with standard deviation of max |<?j'j |) is added to simulate 
measurement error. The conductivity a' G L°°{D') is defined so that a' = cr in 
D'nD and ct' = 1 in £>' \ D. 

Reconstruction is then carried out as above but in the incorrect (or ideal) domain 
D and for the incorrect electrode positions Xj, as if <;-j = <,ij. This simulates the 
effect of slight misplacement of electrodes and error in modeling the domain (while, 
however, still having true relative data). A lower order of M = 12 singular values 
is used. 

Figures 2b-2f show the reconstructions from five different noisy samples. Geom- 
etry error and electrode misplacement are illustrated as described under Figure 2b. 
The effect of these errors, and the extra noise level, in data is visualized in Fig- 
ure 2a, which matches the noisy data corresponding to Figure 2b. This is done with 

the aid of sweep data [15], a restriction of bisweep data where one variable is fixed. 



In Figure 2a, the dashed fine is the sweep data <;a{xi,^]^^{e^^)), £ [0, 27r) of the 
ideal domain D, and it corresponds to the data on a certain vertical (or horizontal) 
line approximately in the middle of Figure lb. The solid line is <;cr'{x[,^~jj}{eJ^)) 
and it illustrates the effect of the geometry error. On the solid line, the values of 
6 corresponding to the actual electrode positions x', are marked and the matching 
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values of c,^' arc the discrete data without extra noise. In addition, the discrete noisy 
data (^[j are plotted at values of 9 corresponding to the ideal electrode positions Xj. 

Approximate inclusion locations are recovered. Also notice the packing of elec- 
trodes near the non-convex corner, which poses an obvious problem for the appli- 
cability of this method in practice. 

5. Discussion and concluding remarks 

The concepts of point electrodes [17] and, in particular, bisweep data [21], 
were generalized from smooth to piecewise 't^^'°' smooth plane domains. This was 
achieved by introducing a generalized relative Neumann-to-Dirichlet map, whose 
properties were studied with the help of conformal maps, Laplace equation with 
measure Neumann boundary values [27], and distributional Neumann problems in 
smooth domains [21] [18]. 

New partial data results for Calderon problem were obtained in Theorem 2.3 and 
Corollary 2.4, based on [21] and the full data result [1]. Instead of assuming a priori 
(interior) smoothness from the conductivity a (as in [25]), boundary homogeneity 
(1.3b) is assumed. In some applications, (1.3b) may be well-founded. For instance, 
if homogeneous medium is screened for hidden defects, it might be reasonable to 
assume boundary homogeneity, but a priori assumption on the smoothness of the 
defects (such as cracks or air bubbles in concrete) could be unrealistic. 

The regularity assumptions for the considered plane domains D is due to Theo- 
rem 3.5. The proof technique based on Theorem A. 2 and weak solutions u in W^'^ 
appears to fail in less regular domains. The other theorems seem to remain valid 
in general Lipschitz domains if the weak conductivity equation is defined appro- 
priately. In this paper, only isotropic (i.e., scalar) conductivities were studied, but 
there seems to be no reason why the results presented here would not have useful 
counterparts in the anisotropic case too. 

It was also demonstrated how point-electrode-based methods could be used in 
numerical EIT, and the notion of bisweep data enables applying unit-disk based 
algorithms in piecewise smooth (polygonal) domains. In this paper, only the factor- 
ization method was considered, but the same approach could be applied to other 
methods as well. It should be noted that there are also other means of applying the 
factorization method in non-smooth domains, for example, [28]. 

The applicability of these point-electrode-based numerical methods depends on 
the availability of relative data [17], or very accurate knowledge of the background 
map Ai , both of which are questionable assumptions in practice. However, theoret- 
ically less-studied concepts such as frequency-difference measurements [19] could be 
used to overcome the problem. This involves analysis of the conductivity equation 
with complex a and such counterparts of the results presented in this paper are left 
for future studies. 
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(a) The finite element mesh used in the 
simulation of bisweep data. 




(c) Unit disk phantom a = a o ^ ^ and 
electrodes xi, . . . , xi28- 




(b) Discrete bisweep data iij, i, j = 1, . . . , 128. 




(d) Reconstruction. The actual boundaries 
of the inclusions are marked with solid lines. 
The different shades of gray correspond to 
different Coo- Electrodes Xi,...,xi28 are 
marked with dots. 



Figure 1. Simulation of and reconstruction from noiseless dis- 
crete bisweep data. 
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(a) Sweep data (solid line, dots (b) Geometry: (?D' and inclusion boundaries 

sampled at Xj, • • • square at x'^), sweep (solid), dD (dashed); electrodes: xi, . . . ,Xn 

data ?ct(^1i-) (dashed), and {xj,^^^) (x). (circles), (square), and Xj, xjj (dots); 

The noisy data equals that of Fig. (b). and reconstruction (gray levels) 




(E) (F) 

Figure 2. (b)-(f): Reconstructions from noisy data. The different 
shades of gray correspond to different Coo (log- spacing). 



22 



OTTO SEISKARI 



Appendix A. Conformal maps in piecewise smooth domains 

Theorem A. 2 shows how any piecewise "^^'^ smooth plane domain can be con- 
formally mapped to the unit disk B through an intermediate (piecewise smooth) 
domain D' so that the maps D' D and D' ^ B are "^^'^ smooth. This is 
based on how piecewise smooth domains transform under fundamental conformal 
maps z !->■ 2", as stated by Lemma A.l. Lemma A. 3 justifies the use of conformal 
transplantation with weak Neumann problems. 

Lemma A.l. Let a S (0,2) and D he a piecewise ^^'^ smooth plane domain such 
that € dD and the mapping z ^ z°' (with some branch cut) is continuous and 
injective in D. 

The mapped domain _D" is piecewise '^^'^ smooth for j3' — min(/3, /3/q;) > 0. 
Moreover, if D has an internal angle of ^ ai and a = ^, then is "if ^■'^ smooth 
near the origin. 

The above can be proven easily using a natural parametrization of the boundary 
dD. 

Theorem A. 2. Let D be a piecewise ^S"^ " smooth plane dom,ain and <I> : D — > 
B a conformal map. There exists a piecewise 'Tf^'°' smooth plane domain D' and 
conformal maps : D' ^ B, ^2 ■ D' ^ D such that $1 e "^^'"'(Z) ) for some 
a' > 0, $2 e '^°°(d') and $ = $1 o 

Proof. The domain D' can be constructed as follows: pick a corner x G dD with 
a reflex internal angle 7 > tt (if any) and some point x' outside D such that the 
line segment between x and x' does not intersect D. The conformal map z i-)- 
l/{z — x') G '^°°{D) yields a piecewise smooth domain Di with an internal 
angle 7 at 5 = l/(.x — .x'), and there exists an open ray from x to infinity that does 
not intersect Di. Using this ray as a branch cut, one may apply the transformation 
z ^ {z — .f)^/^, which yields a domain that is piecewise "^^'^ smooth and has 
one reflex angle less than D. This process may be repeated until no internal angle 
is reflex. The resulting domain is D' , and the chained (inverse) transformation $2 
is '^°°{D') by explicit construction. 

A conformal map ^ : D' ^ B can be constructed as a composition $3 o $4 =: 
$3 ° $m o • • • o ^1, where : D'^_-^ D'j, 

$,(z) = (</>,(z))-/''^- , e ^°°(5^), .7 = 1, . . . , m, 

straighten the angles 7^ < % oi D' =: D'q. Thus, for some a' > 0, $4 G "^^'"'(D'), 
and the map $3 G "^^'^ {D'^) transforms the resulting smooth domain to the 

unit disk B (cf. [35, §§ 3.3 & 3.4]). 

Finally, any two conformal maps $ : D — B are related by ^ = Mo $, where 
M G '^°°{B) is a Mobius transformation, which shows that $ = M o $ o $2 ^ ='■ 
$1 o $2 ^ is as claimed. □ 

Lemma A.3. Let f,ge Hl^^{D) be such that Vf-VgG L^{D). Define / = / o 
g = g o ^ and (J = (7 o where ^ : D ^ D is conformal, D, D Lipschitz, and 
d G L°°{b). Then Vf-Vg€ L^(D) and 

/ aVf - Vgdx = / aVf-Vgdx. 

J D J D 
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Proof. Let = {x e D : dist{x,dD) > e}. Clearly <I>\d, € -^""(D,) for all e > 
and thus 

[ \(7Vf-Vg\dx= [ |(cto$)V(/o$). V(5o$)|da; 

= [ |(ao$)(j^j|(V/)o$).((V5)o$)|dx 



(A^l) 



= / 

-j 



a.(r>e) 



(det Ja, o $-^)ctV/- V5 
d-Vf-Vg\dx, 



I det J$-i| da; 



where J<^ : D ^ M?^"^ denotes the Jacobian of $ (interpreted as a mapping D c 
R2 ^ $ conforniaL J<i,jJ = = (dct J<i,)/ = (dot J$-i o $)-iJ. 

The monotone convergence theorem yields that that ||crV/ • yg\\L^{D) = ll'S'V/ • 
Similarly to (A.l), it holds that Jj^ aVf ■ Vgdx = J^f^jj ) ctV/ • Vgdx 
and the claim follows from Lebesgue's dominated convergence theorem. □ 
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